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Abstract. The physics of the attractive one-dimensional Bose gas (Lieb-Liniger 
model) is investigated with techniques based on the integrability of the system. 
Combining a knowledge of particle quasi-momenta to exponential precision in the 
system size with determinant representations of matrix elements of local operators 
coming from the Algebraic Bethe Ansatz, we obtain rather general analytical 
results for the zero-temperature dynamical correlation functions of the density 
and field operators. Our results thus provide quantitative predictions for possible 
future experiments in atomic gases or optical waveguides. 
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1. Introduction 

The experimental realization of trapped quasi-one-dimensional atomic gases in optical 
lattices during the last few years [1, 2, 3, 4, 5, 6, 7, 8] has provided a new impetus for 
the study of the effects of strong correlations on the physical properties of fundamental 
quantum-mechanical systems of interacting particles. One very interesting aspect of 
this still emerging field, in contrast to traditional condensed matter physics, is that 
more-or-less ideal realizations of famous toy models like the interacting Bose gas in 
one dimension (the Lieb-Liniger model [9]) can be constructed and investigated. It is 
moreover usually possible to tune the parameters to arbitrary values, in contrast to 
the fixed values in a solid-state crystal, and to obtain systems in which disorder and 
impurities are essentially absent. 

Experimental measurements can typically be related to relatively simple 
theoretical objects like expectation values and correlation functions of local fields. 
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Ideally, one would like to obtain direct quantitative comparisons between these, 
providing an extremely nontrivial check on our understanding of the physics involved. 
From the theory side, calculating the response functions of a strongly correlated system 
is more often than not a difficult, if not insurmountable task. In one dimensional 
quantum mechanics, the situation is not as hopeless as in higher dimensions in view 
of the existence of many nonperturbative methods [10]. While features such as 
asymptotic correlation functions can be obtained from bosonization [11, 12], it remains 
however extremely rare to be able to provide general analytical results valid for all 
ranges of momentum and frequency. 

Our objective in this paper is to provide an extensive treatment of a situation 
where much analytical progress is possible, namely that of the Licb-Linigcr model 
for a one-dimensional Bose gas in the attractive regime. Like the repulsive case, 
the attractive one-dimensional gas is solvable with a Bethe Ansatz, so we can hope 
for full nonperturbative solutions. On the other hand, although the Bethe Ansatz 
offers a good basis for providing an 'exact' solution to some models of interest in 
one-dimensional quantum physics, it is not necessarily an easy framework to use 
to compute dynamical correlation fimctions or other objects going beyond simple 
equilibrium thermodynamics. In the case of the Lieb-Liniger model, while it is now 
possible to obtain accurate numerical results for the dynamical structure factor [13] 
and for the one-body function [14] for the repulsive Lieb-Liniger model by mixing 
the Algebraic Bethe Ansatz with intensive numerical work, obtaining a full analytical 
solution for these objects in the thermodynamic limit remains a severe challenge [15]. 
Here, we demonstrate that the attractive gas is in fact fully tractable analytically. A 
summary of our results has appeared recently in [16]. 

The physics of the attractive Lieb-Liniger model, less conventional than in the 
repulsive regime, has received attention in the past. Most importantly, the ground 
state energy and wavefunction were investigated in [17, 18], and a more general 
treatment including excitations was offered in [19]. The purely one-dimensional case 
also serves as a basis for the investigation of effects beyond the Gross-Pitacvskii mean- 
field in higher dimension [20]. Overall, the scattering and diffusion of solitons has 
also been studied [21, 22], and mean-field approaches based on the Gross-Pitaevskii 
equation have been developed [23, 24]. The case of attractive interactions is in principle 
accessible experimentally in the context of bright solitons in quasi ID traps [25, 26, 27] 
since the effective one-dimensional coupling constant [28, 29, 30] can effectively be 
tuned via Fcshbach resonance or transverse confinement to essentially any positive or 
negative value. Such sort of solitonic physics can also be realized using photons in 
optical fibres [31]. 

The paper is organized as follows. We start by a reminder of the model, its 
solution in terms of a Bethe Ansatz, and of the important correlation functions to 
consider in Section 2. After describing eigenstates in more detail in Section 3, we 
give the important form factors in Section 4, which are then used in Section 5 and 6 
to obtain the dynamical structure factor and one-body function at zero temperature. 
Most of the more tedious calculations have been relegated to Appendix A- Appendix 
E. We end by offering some conclusions and perspectives in Section 7. 
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2. Setup 

The Hamiltonian of the Lieb-Liniger model is given by 

where {i,j) represents the sum over all pairs, c is the interaction parameter and m the 
mass of the particles (atoms). In terms of experimental parameters the ID coupling 
constant is c = —h^/mauj, where auj is the effective ID scattering length that can be 
tuned via Feshbach resonance or transverse confinement [28, 29, 30]. For definiteness, 
we consider a system of length L (i.e. < xt < L) with periodic boundary conditions. 
From now on we fix ?i = 2m =1. In second quantization, this is nothing but the 
nonlinear Schrodinger theory for a canonical Bose field, 

H = [ dx{da^¥{x)da^^{x)+c^\x)^''{x)^{x)<b{x)} . (2) 
Jo 

Wc will be primarily interested in obtaining detailed results on dynamical correlation 
functions of local operators of physical significance for the attractive gas (1) with 
c < 0. We will first consider the density operator, which is defined as 

^ 1 

p{x) =J2^ix- xj) = j^e^'^p, = ^\x)^ix). (3) 
j=i k 

More precisely, we will obtain its ground-state correlation function| 

S''{x,t) = {pix,t)p{0,0)), (4) 

whose Fourier transform S''{k,oj) is known as the dynamical structure factor (DSF). 
In a similar way, we will also treat the case of the one-body dynamical correlation 
function of the field operator ^{x,t), 

S^{x,t) = {^Hx,tMO,0))- (5) 
The dynamical structure factor (4) can at least in principle be obtained 
experimentally by Fourier sampling of time-of-flight images [32] or Bragg spectroscopy 
[33, 34]. The one-body function (5) is obtainable through interference experiments 
[35, 36] or Raman spectroscopy [37, 38, 39, 40, 41], and its static version is the 
momentum distribution function which is directly accessed using ballistic expansion 
[42, 4]. 

By introducing a resolution of unity, these correlation functions can be written 

in a Lehmann representation involving a sum over intermediate states. As a basis for 
the Fock space, we use the eigenfunctions of the Hamiltonian as given by the Bethe 
Ansatz, which are combinations of plane waves characterized by spectral parameters 
(rapidities) {A}: 

*jv(a;i, ...,a;jv|Ai, Aiv) = Y\ sgn{xj - Xk) x 

N>j>k>l 

Pn 

t Every time we write an operator with an explicit time dependence we obviously are working in 
the Heisenberg picture. Oppositely operators without time dependence (as in Eq. (2)) are in the 
Schrodinger picture 
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where (j){X) = 2arctan(A/c). These wavefunctions are symmetric in the position 
coordinates. The rapidities of a given cigcnstatc obey a Pauli-Ukc exclusion principle, 
namely the Bethe wavefunction vanishes if two rapidities coincide. Quantization is 
achieved by enforcing periodic boundary conditions. It will become clear later on that 
our results are in fact insensitive to the boundary conditions used. 

Upon Fourier transformation 0{k,u)) = dx J^^dte^'''^*~''^^0{x,t), the 
correlators we are interested in can be written as a sum over intermediate states 
labeled by jj., 

S'^{k,u;) = 2^L^ 6{u; ^ E, + E,o), (7) 



where A° denote the ground state rapidities, and the form factor (FF) = 
(/x|O(0, 0)|A) = (/u|Oif^_if^|A)/L depends on the operator O and on the state n {Ok 
is the Fourier transform of the operator O at momentum K). Throughout this paper, 
denotes the norm of state fj.. 



3. Eigenstates 

Let us begin by discussing the eigenstates, and by first recalling some standard features 
of the Bethe Ansatz (6) for the model at hand, and how they adapt to the attractive 
case. 

As can be seen from the Schrodinger equation, given a set of rapidities, the energy 
and momentum of the corresponding state are 

N N 

Quantization of the theory is achieved by imposing the periodicity conditions 

*(a;i,a;2, ...,a;jv) = *(a;2, ...,a;jv,a;i + L) (9) 
which lead to the Bethe equations 



e 



or in logarithmic form 

Aji + ^.^(A,-Afc) = 27r7j, j = l,...,N. (11) 

k 

The quantum numbers Ij are half-odd integers if N is even, and integers if N is odd. 

In the repulsive case, given a proper set of quantum numbers {/}, the solution 
of the Bethe equations for the set of rapidities {A} exists and is unique [43] due to 
the convexity of the Yang- Yang action associated with (11). Furthermore all these 
solutions have real rapidities A^. For the attractive case, the situation is completely 
different. We will define c = — c > as the interaction parameter; let us rewrite our 
Bethe equations as 

ilA„-A, + zc- 
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Figure 1. Two string states of 
a gas o{ N = 7 atoms. Black: 
The ground state consists of a single 
string centered at fc = with all the 
N particle aligned on the imaginary 
axis. Red: An excited states with 
4 strings of length j = 1,2,3 and 
Ni = 2,N2 = 1, Ns = 1, Nj^3 = 0. 



Consider now a complex rapidity Aq, = A + irj. The Bethe equation for this rapidity is 

J- J- A, - A/3 + zc 

We consider finite N and L oo. If 77 > 0, we have e^"^^ ^ on the left-hand side. 
Looking at the finite product on the right-hand side, we conclude that there must thus 
be a rapidity \a' such that Aq.' = Xa — ic + 0(e~''^). On the other hand, if 77 < 0, 
we have e^''^ cx3 on the left-hand side, and there must thus be a rapidity Aq' such 
that Xa' = A„ +ic-|-0(e-l''l^). 

In general, the atoms will thus form bound states in the form of strings in the 
complex plane. For a given number of atoms N, we can construct eigenstates with 
fixed string content by partitioning N into Nj strings of length j, denoting the total 
number of strings as Ng. We clearly have 

iV-^jW,, N, = J2^r (14) 

j 3 
Specifically, we will parametrize the string rapidities as 

Ar = Ai+7^(j + l-2a) + z^r, (15) 

with exponentially small deviations 5 ^ g-i'^'^t)^ as usually done in these cases (see 
e.g. [44, 45]). In our string notation, the index a = 1, . . . ,j labels rapidities within the 
string, and a = 1, . . . ,Nj labels strings of a given length. A pictorial representation 
of two string states is given in Fig. 1. We stress that perfect strings (i.e. with all the 
6i = 0) are exact eigenstates in the limit L ^ 00 for arbitrary N. It is then natural 
to consider the limit L — s- 00 at fixed N. This is different from what done in the 
repulsive case where the limit iV, L ^ 00 at fixed density N/L is performed. Here, the 
N particles remain strongly correlated and bound to one another even when L 00, 
as we will explicitely see later on (see in particular the section about the norm of a 
string state 3.3). 

The string deviations 5^'" are effectively studied in [46, 47], but depend sensitively 
on the particular boundary conditions used. We will however not need the explicit 
values of these deviations: we will only need the fact that they are small, and treat 
the resulting limit carefully in the norm and form factor expressions. For the moment, 
however, they can be dropped. 

These bound states should be viewed as individual particles of mass j, with 
momentum and energy of the string centered on A^ given by 

Eu,c.) = jiK)' - ^Jif - 1), = jK- (16) 

Such strings are known but not commonly discussed in the literature on the Bose gas 
[44, 19], as they do not appear in the repulsive case. However, their equivalents exist 
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in integrable spin chains, where they have been extensively studied. The 'technology' 
to treat them, at least on the level of cigcnstatcs, is thus standard. 

The fact that all deviations are exponentially small means that we can obtain a 
much simpler set of equations involving only the string centers A^. The derivation of 
these equations in presented in Appendix A, the final result being 

jXiL - ^ <l>,,(Ai - A^) = 2nli (17) 

(fc,/3) 

with half-odd integer (integer) if Nj is even (odd), and where the scattering phase 
shifts are 

*,fe(A) = (1 - S,k)4'ij_k\W + 20|,_fe|+2(A) + ... + 2(/.,+fe_2(A) + (f'j+kW, (18) 
with 0j(A) = 2 arctan(2A/jc). These strings are stable particles under scattering with 
one another, and are therefore soliton-like objects. One point worth emphasizing is 
that the scattering phase shifts (18) arc simply those of breathers in the classical limit 
/3 ^ of the sine-Gordon model after a trivial reparametrization of the rapidity. The 
sine-Gordon soliton's mass is in this limit much higher than that of the breathers, 
which have an evenly-spaced rest mass. The string of order j in the attractive 
Lieb-Liniger model is consequently not so much a 'soliton', but more accurately a 
nonrelativistic sine-Gordon breather of order j. 

3.1. The ground state 

The lowest energy state will be obtained by forming a bound state of all N particles 
centered on zero [17] (see Fig. 1), namely by choosing 

X^^" = i^{N + l-2a) + 0{8). (19) 
The corresponding energy is 

-2 ^ -2 

Egs = = + 1 - = -^^(^' - !)• (20) 

a a=l 

Note that this goes like , imlikc the repulsive case where the energy is not extensive. 
Although our results are not limited to this case, we will often consider the limit of 
a large number of particles N ^ 1, with weak interactions c such that the parameter 
g = cN remains finite. In this case, the ground state energy per particle also remains 
finite, Egs = —5^/12. We will find some similarly simplified limiting values for the 
correlations functions. It is however important to note that this is not a conventional 
thermodynamic limit with finite energy density as in the repulsive case. 

3.2. Excited states 

Excitations above the ground state are then obtained by either simply giving 
momentum to the ground state N string, or more elaborately by partitioning it into 
smaller strings to which individual momenta can be given. We consider here only the 
N atom sector relevant for the dynamical structure factor. The case of states with 
A'' — 1 (or less) atoms trivially follows. 

We will label the string content of eigenstates by column-separated entries 
specifying the length and number of each different string type. For example N—M : M 
will be a state with a A^ - M string and a M string, and N - Mi - 2M2 : Mi : (M2)2 
a state with an A^ — Mi — 2M2 string, an Mi string and two M2 strings; the red state 
in Figure 1 is a 3 : 2 : I2 state. 
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3.2.1. Single-particle states Single-particle excited states will be obtained by giving 
finite momentum to the ground state A''-string, 



^^.« = ^ + i|(7V+l-2a) + 0(^). (21) 



Such states have energy above the ground state given by 

5^ - Egs = Nfi^ = kl 



CON (m) = E^,- Egs = Nfi^ = kl/N (22) 



where /c^ = N fj, is the total state momentum. For these states, there is only one Bethe 
equation for the string center fi, namely /i = 21^1 /NL with I an integer, so that the 
momentum is quantized as for a free wave, = 2-^1 jL. In the limit of large iV, this 
energy band becomes flat and quasi-degenerate with the ground state. 

3.2.2. Two-particle states These are obtained by splitting up the groimd-state N 
string in two pieces. In general, consider having an TV — M and an M string; 

^N-M,a = + i£(7v - M + 1 - 2a) + 0((5), a = 1, N - M, 

H^^"- = fXM + i^{M + l-2a) + 0{5), a = l,...,M. (23) 

The energy of this state above the ground state is given by 

iON-M:M{lJ's,I^M) = l^%-M:M + {N - M)fll + Mj^lj, (24) 

where we have defined the rest energy 

to%_M..M^^NM{N-M). (25) 

The total momentum is the sum of the two string momenta, 

k = ks + kM = {N- M)ij.s + Mum, (26) 
so we can write the energy as 

UJN-M:M{ks,kM) = i^N-M:M + ^ + -JJ ■ (27) 

Similarly to the single-particle case, the Bethe equations are here very simple, namely 

{N - M)^isL - ^N-M,Mil^s - fJ'M) = 27r/s, 

MiimL + <^n-m,m{i^s - I^m) = 27r/M, (28) 
with 7s, Im integers. In the limit of large L, we can thus ignore the scattering phase 

shift, and take /i^ and as free parameters. The total momentum k of the state 
can take on any value 2tiI / L, but the energy is bounded from below by 

I 

^n-^m-.mQ^) = ^N-M:M + j;;^- (29) 

Given external frequency w and momentum k parameters, there are two solutions to 
the dynamical constraints, namely 



M 



_i_ k 
MM(fc,w) = ;^ ± 



N{N - M) 
N-M' 



1/2 

[w-w5v-M:M(fc)]'/^ 



1/2 

j^j^j [^-^'N-M:Mik)Y^^- (30) 

Therefore, in the large L limit, these states for a two-fold degenerate continuum 
beginning at the lower threshold i^N-M-mi^) ^^'^'i extending to arbitrarily high energy, 
N-M:M continuum: lj^j^ _M.M{k) < lo < oo. (31) 

For finite L, this is of course not strictly a continuum: only discrete energy levels u> 
then exist, as determined from the Bethe equations. 
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3.3. Norm of a general state, general remarks 

Before discussing form factors and correlation functions, we need to complete our 

characterization of individual cigcnstatcs by providing a formula for their norm. Given 
a set of rapidities solution to the Bethe equations, the norm of the corresponding Bethe 
eigenstate is given by the Gaudin-Korepin formula [48, 49], in this case 

ll*iv({A})|p = |cr n %^Detjve (32) 

j>k ok 

where Q is the Gaudin matrix whose entries are given by 

N 

5jfc({A}jv) = 5jk 



K{Xj,Xk) (33) 



where the kernel is 



To cover the case of string states in the attractive regime, these formulas must be 
adapted in view of its string deviation 6^^°- / S^^"-^^ indeterminacy. We do this in detail 
in Appendix B, where we obtain the final result 

ii*iv({A^j)ip=(ic-)^'nj"'^^ n pM->^''0) (35) 

3 (fc,/3)>(j,c<) 

where 

^^■'=W-A2 + (|(,-fc))2- (36) 
The norm of a single N string is thus simply 

||A||2 = cLiV2 (37) 

and that of an iV — M : M state with rapidities /Ltg, jim respectively for the N — M 
and M strings, 

\\{l^s,m}\? = c'L^N - MfM^^,^ ^ + ^ ^ (38) 

As we already mentioned, unlike what happens in the repulsive case, where the 
low-density limit N finite, L ^ oo is trivial, the physics of the attractive case remains 
nontrivial. This is not a low-density limit, as the atoms tend to clump together in 
wavepackets of finite extent. The ground state wavefunction is thus really a localized 
wavepacket spread out uniformly over the whole system, as can be seen from our choice 
of normalization: this depends on L generally as i^", thus counting the number 
of strings present, reflecting the fact that the strings are essentially independent, 
almost free particles. In the presence of a weak harmonic confining trap in an 
experiment, the ground state would thus be obtained by a convolution of the N string 
wavepacket wavefunction with a ground state harmonic oscillator part. Our results 
could straightforwardly be adapted to cover this case. 
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4. Form factors 



Since the Bethe wavefunctions are rather involved, comprising a factoriaUy large 
number of separate free wave terms, it is a priori extremely difficult to calculate matrix 
elements of local operators in this basis. However, the calculation can be performed 
for a number of integrable models (including the Bose gas) using the technology of 
the Algebraic Bethe Ansatz, which is sufficiently powerful to provide a handle on 
this combinatorial- like problem. In fact, the ABA provides these matrix elements for 
the density and field operators in terms of determinants of matrices whose entries are 
simple analytical functions of the rapidities of both eigenstates involved [50, 51, 52, 14]. 
These formulas are also applicable, at least in principle, to the attractive Bose gas, 
provided some care is taken with string deviations. The purpose of this section is to 
perform this calculation. 

For the repulsive gas. the form factors have to be evaluated numerically, since the 
solutions to the Bethe equations cannot be found analytically. For the attractive case, 
however, we can actually calculate these determinants analytically to leading order 
in the exponentially small string deviations, since we know the ground state to this 
precision without having to solve any Bethe equation. 

The form factors of the field and density operators between the ground state and 
an arbitrary excited state are computed in Appendix C and Appendix D, and are 
given by 



p2 

c 



c^N\{N ■ 



where we have defined the function 



Hm{x) 



-l)\X{HMc), 

1)! n HM^)^ 

(3,01) 



ix) 



T{^+ix) 

For later convenience, the asymptotic properties of H as N ^ M are given by 

-^2 n-M 



(39) 



(40) 



(41) 



Hm{x) 



N-M 



{X) 



Fm{x) 



(l + 0(M/iV)), 



[N^+x^]""^ {1 + 0{M/N)), 



F2n{x) 

where we have defined the function 

Fm{x) = \T{M/2 + ix 
For X <^ N, we have further that Hn-m 



(42) 



(43) 



{x) ^ N^FM{x)/r\N) + 



5. Density correlation function 

Let us now apply the results previously obtained for states and form factors to the 
calculation of the dynamical structure factor. We first present the sum rules by which 
our results can be checked, and then treat each important class of excited states 
separately. We also consider how the equations obtained simplify in the limit of many 
weakly interacting atoms. 
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5.1. Sum rules 

The first sum rule is for the second density moment (p^), which is given by the 
Hellmann-Feynman theorem as 

and we must therefore have 

lE/, -^s<'{k,^) = {^p{mp{m) = ^-^ — -■ (45) 

A: 







Another sum rule of practical interest is obtained from the first frequency moment at 
fixed momentum (/-sumrule): 

r(fc)= '^u:S\kM = ^k\ (46) 

We will saturate these sum rules in what follows by adding up contributions to the 
structure factor coming from the lowest-lying bands of states. 

5.2. N contribution (one-particle) 

When the rapidities {/i} correspond to an A'^-string, we have from Eq. (39) that 

1^, ^ r(zf )r(-zf )r(jv + i)r(7V) p'^ _ Nr^iN)p^^/c 

' r(iv + zf)r(7v-*f) c Ua^,[(p/cr + {a-m ^ ' 

and k = P^ = N^. The norm of this state is ||{m}|P = cLN^, so 



S'^{k,u,) = 2n— }J'l, S{u; -). (49) 



This contribution is thus a single peak centered on the N string dispersion relation. 
At large A'' with g = cN fixed, this becomes 

N^Jnk/g)^ e 

' 2 

L sinh nk/g N ' 

The density moment sum rule contribution from these intermediate states is for 
large N: 

r-5^(fca;)-^ T-^^^^-— (50) 
i^^io 27r^^^^''"^^ L J_^2nsmh'ink/g)- 6L ^^^^ 

so these states completely saturate this sum rule for large A''. 
The /-sumrule contribution is simply written as 

U2 



which for large N and fixed g becomes 



nk/g 



< ^e. (52) 



smh-Kk/ g 

These states therefore do not completely saturate the /-sumrule, so we consider higher 
excited states in order to find the missing part. 

Note that if we instead considered ^ oo at fixed c, the /-sumrule would have 
been saturated by these states, i.e. the dynamics would have been less interesting. 
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5.3. N — M : M contribution (two-particle) 

The next simplest excited states to consider are the two-particle N — M : M states 
with M = 1, 2, .... The rapidities of these arc parametrized using (23), and the energy 
and momentum with respect to tlie ground state are respectively given by (24) and 
(26). Given an external frequency w and momentum k within the continuum (31), 
there exist two solutions to the dynamical constraints for the rapidities iisi^-M given 
by equation (30). 

From (39), we find that the density operator form factors relevant to this class of 
states are 

p2 

|E;;^o| = Nr\N)^HN_M{f^s/c)HM{m/c). (53) 

Their contribution to the dynamical structure factor is thus 

IS'' P 

S'N-M:Mik,U;) = ^ n , p^, 1 12 '^^.^/('^ " ^/^ + ^A")- (54) 

Let us explicitly rewrite the right-hand side as a function of the variables k, w. The 
energy delta function can be rewritten using (30) as a function of /Xgj 

M Y'^ e(w-a;^_M:M(fc)) 



N{N -M) 



Y^8{iXs-ifAkM)-{^^) 

o-=± 



For large L, the summation over /ig can be replaced by an integral, — > L 
[iM is then completely fixed, finally yielding 

. e(u;-u;5v-M:M(fc)) V\N)k^ 



As everywhere, this expression is correct to leading order in the exponentially small 
string deviations for any value of N and M. Each choice of M gives a distinct 
continuum. 

In the limit N ^ M 1, this becomes 



@{UJ - UJ^^_M,M{k)) 

X 

J^m(m?/c) 



'^N-M:M\'^^^) — ^ , Nnl/2 



■ 2iV2M+3M3/2c5L [1/4 + (^^/c7V)2: 



2M ■ 



(57) 



The hierarchy is such that higher M parts are suppressed by increasing powers of iV§. 
The leading M = 1 part can be further simplified to 

*jv-i:i(fe>'^)- i25[w-5V4]W ^^^^ 
with Fi{x) = 7r/cosh(7ra;) and x± = ([a; — 5^/4]i ± k)/g. 

§ If M is of the order of AT it is straightforward to show that the contributions are exponentially 
suppressed for large N 
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SP ,^, (k,0))/(N/L),g = 1 




Figure 2. Left: contribution to tlie DSF coming from A'^ — 1 : 1 (two-particle) 
states, in the large TV limit for g = 1. The square-root singularity at the lower 
threshold is accompanied by a maximum around k = [u) — 2 . Right: Fixed 

momentum cuts of Fig. 2, showing the weakening of the square-root singularity at 
higher momentum, and the displacement of the maximum towards higher energies. 



The contribution of these state to the / sum rule is written as the integral 
(introducing gC = ^Juj ~ g'^ /A) 

' " L y_ 27r e + 1/4 cosh^ 7r(C - k/g) ' ^''^ 

In the upper half-plane, the first piece of the integrand has a single pole a,t ( — i/2 
and the second double poles at C = k/g + i{n + 1/2), n = 0, 1, .... Writing out the 
residues gives a series which can be explicitly resummed, leading to the final result 

.X NkW (nk/g)^ \ , , 



i \ sinh (jikj g') ) 

These states therefore completely saturate the remaining part of the /-sumrule after 
the iV string contributions (52) have been taken into account. Higher frequency 
moments further suppress the S"^ part; in general, any frequency moment can be 
accurately obtained by approximating the dynamical structure factor with 5''(fc, iS) — 
^^(fc,^) + S''^_^..y(k,uj). 

The two-particle (iV — 1:1) part of the DSF is plotted in the left panel of 
Figure 2 for large N , with g = \. At the lower threshold, the DSF diverges as 
(oj — 9^/4)^^/^. For Lj > g'^ /A it is a monotonous decreasing function of lo as long as 
k/g < Xc = 1.0565 . . ., whereas for k/g > Xc it shows a broad peak, whose position 
grows like fc^ for large k (its amplitude decreases like fc~^). Away from this peak 
it decays exponentially. The right panel of Figure 2 provides fixed momentum cuts 
showing these features in more detail. 

5.4- Static structure factor 

As a byproduct of the previous results we can simply obtain the static structure factor 

/•OO 7 

S'ik)^ ^S^ik^u) (61) 
Jo 27r 




2 3 45 012 3 45 

k/g k/g 

Figure 3. The static structure factor. Left: Two-particle contribution 
S'^_j.j(fc)/(Af/L). Right: One- plus two-particle contribution (S^(fc) -|- 
5^_j.j(fc))/(Af2/L) for N = 10,100,1000,100000 showing the crossover at 
k/g ~ InAf. 



from our results. The contribution of one particle states is given by the integral of 
Eq. (49) over uj, yielding 



^ L sinh^(7rfc/(7) 



S''N{k) = — ^ . (62) 



We can also evaluate the contribution of the N —1 : 1 states 
Once again, the residues series can be explicitly resummed, giving 



smh^TT- V 9 



(64) 



where ^2{z) is the polygamma function (second derivative of the logarithm of the 
Gamma function). A plot of 5^_-^.-^(fc) is shown in the left panel of Fig. 3. It shows 
a broad peak at k/g ^ 1.2. This is similar to what observed numerically in the 
super-Tonks-Girardeau gaslike regime in Ref. [53]. Although any eventual connection 
between the two results is far from being clear, it is tempting to believe that the 
quench "experiment" used to realize the super-Tonks-Girardeau gas has selected a 
two-particle state, in analogy to what claimed in [54]. 

The one-particle contribution S'^(fc) for large N grows like A^^, and so it will 
always dominate against the two-particle one (growing like N) for any k at large 
enough A^. However for large fc/g, S'/^{k) decays exponentially, whereas '5'^_x i(^) 
goes to the constant value N/L. Consequently, a crossover between the two regimes 
is expected to take place at relatively low momentum even for very large A^. In fact, 
imposing >5'^(fc) ~ N/L we have kc/g ^ (In A^)/27r. This crossover is explicitely shown 
in the right panel of Fig. 3. 

6. One-body function 

For the dynamical one-body function S^{k,uj), the relevant intermediate states are 
those made of A^ — 1 atoms. When calculating the energy difference of a given 
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excited state with the ground state, we will thus have to take the chemical potential 
E% - E%_^ ^ -c^N{N - l)/4 into account. 

The relevant sum rule here is very simple, namely that we should recover the 
particle density by integrating over energy and summing over momenta, 



6.1. N-1 contribution (one-particle) 

In this case, we write the excited state as an iV — 1 string, 

/x^-i.«=;x + i|(7V-2a), a = 1, TV - 1. (66) 

The total momentum of this state is P = (A'^ — and its energy above the 
ground state (taking the chemical potential into account) is (iV — l)/x^. Its norm 
is = cL{N — 1)^, and its field operator form factor with the ground state is 

obtained from (40) as 

l^U - „ . . (97) 

n.."i iwc)^+(«-i/2)'i 

The one-body function contribution is thus 

In the large A'' limit at g constant, this becomes 

27r^A^ 5{oj - k'^/N) 
gL cosh^ {wk/g) 
Integrating this over frequency and momenta, we find that 

_f ^5*_i(fc,a;) = ^ (70) 

so these states completely saturate the integrated intensity sum rule in this limit. 
Since we can also be potentially interested in higher frequency moments, which would 
then not be completely saturated by this function, we consider higher excited states, 
the leading ones being the two-particle ones. 



6.2. N-1-M:M contribution (two-particle) 

In complete parallel to the N — M : M states considered for the structure factor, wc 
here consider the — 1 atom states with defined by having an — M — 1 and an M 
string: 

^N-M-i,a = + ^ c - M - 2a) + 0((5), o = 1, A^ - M - 1, 

M^'° = MM + i|(M+l-2a) + 0((5), a = l,...,M. (71) 

The energy of this state above the ground state (again taking the chemical potential 
into account) is given by 

WjV-M:m(Ms, Mm) = W^_M-l:M + {N - M - l)jUg + Mll\j, (72) 
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where the rest energy is 

uJ%^M-i:M = jiN-miN-M-l). (73) 
The total momentum is the sum of tlic two string momenta, 

k = ks + ku = {N -M - + MfiM, (74) 
so we can write the energy as 

UJN-M-V.M{ks,kM) = i^N-M-l:M + _ M - 1 ^"^^^ 

The Bethe equations are here 

{N - M - l)fisL - $7v-m-i,m(Ms - Mm) = 2irls, 

M/imL + ^n-m-i,m{i^s - I^m) = ^ttIm, (76) 

with Is, Im integers. In the Umit of large L, we again ignore the scattering phase shift, 
and take iig and /zm as free parameters. The energy is again bounded from below by 
^N-M-i:Mi^) — '^N-M-i-.M + k^/N. The dynamical constraints give in this case 

-,1/2 



M 



k 

HM{k,Uj) = j^—^± 



N-1 ' [{N - 1){N - M - 1) 
N-M-1 



[oj - ojf^_]^_i,j^f{k)] 



1/2 

[OJ - OJ^N-M-l-.Mik)]'/^- (77) 



(iV- 1)M 

Therefore, in the large L limit, these states again form a two-fold degenerate continuum 
beginning at the lower threshold '^jv-m-1:m(^) ^^'^ extending to arbitrarily high 
energy, 

N-M-1:M continuum: "^at-m-IimC^) < < oo. (78) 

The calculation of the contribution of theses states to the one-body function 
follows the same logic as that used for the N — M : M states for the structure factor. 
We find 

C,* tl. \- Q{<^-<^N-M-l:Mik)) 

'^N-M-1:M\'^J^) 



1/2 



2cL{N - l)i/2[(7V - M - l)M]3/2 ^ Fjv-m-i,m(m? - Mm) 
At large N, this can be simplified to 

S.-M-.Mik,.)= 1/2 2,2^^3/2 1. [^/4^(^.„/^).j2M- (^0) 

The leading M = 1 term then becomes in this limit 

with again x± = n{y/uj — 5^/4 ± k)/g. Here, ^2(0;) = 7ra;/sinh7ra;. This contribution 
is plotted for g = 1 in Figure 4. 

Note that the relative contribution of these states to the sumrule (65) goes like 
0{N^)/L and so it is suppressed by a factor A'' compared to the one-particle state. 
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Figure 4. Left: contribution to the logarithm of the one-body function coming 
from N — 2 : 1 (two-particle) states, in the large A'^ limit for g = 1. Right: Fixed 
momentum cuts of the same, showing the same qualitative features as for the 
DSF, namely the weakening of the square-root singularity at higher momentum, 
the appearance of a local maximum as a function of frequency for larger values 
of momentum, and the displacement of this maximum towards higher energies at 
higher momentum. 



6. 3. The momentum distribution function 

As for the structure factor, the static hmit, known as momentum distribution function, 

S^ik)^ / pS^ik^co), (82) 
Jo 

is obtained easily integrating the previous results. 
The one-particle contribution is 

Sl-^{k)^ r^SZ^,ik,.u) = ^—^f—- (83) 
Jo '^^ 9^ cosh [nk/g) 

that again decays exponentially for large k. 

We can also evaluate the contribution of the iV — 2 : 1 states to the static one-body 
function for which we get the integral 

dC 1 iC-k/g^ 
gL 2tt (C2 + 1/4)2 gi^i^2 _ j^^g^ ■ 

This integral can be done summing over all the residues. The final result involves 
polygamma functions and it is not really illuminating. For this reason we only plot it 
in Fig. 5. 

Let instead study in detail the behavior for large k. In this case we have 

S.^.,ik » ,) . --^ = (85) 

This fc~^ behaviour is the same as for the repulsive regime [56, 57]. In fact this is not 
a coincidence: The derivation of Ref. [56, 57] relating the large k behavior of S'*(fc) to 
the second density moment (p^) makes use only of analytical properties of the Bethe 
wave function that hold both in attractive and repulsive regimes. Adapting the result 
of Ref. [56, 57] to our normalization we get 

S^ik^^) = {p')^, (86) 
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Figure 5. The two-particle 
contribution to the mo- 
mentum distribution func- 
tion gLS'^_2.i{k). Inset: 
Lk'^S^_2.i(k)/ to sliow the 
k~'^ tail with prefactor 1/6 
(straight line). 



and inserting the value of (p^) obtained via the HeUmann-Feynman theorem in Eq. 
(44) reproduces exactly our last result (85). This is not a trivial result because it is 
saying that the large k behaviour is completely determined by the two-particle states. 

7. Conclusions and perspective 

We presented a detailed calculation of the zero-temperature dynamical correlation 
functions of the attractive Lieb-Liniger model combining the determinant representa- 
tion of the matrix elements with the exact (with exponential precision in the system 
size) knowledge of the Bethe states. All the sum rules we have studied are essentially 
saturated by the very simplest classes of excitations: the one-particle states obtained 
by giving finite momentum to the ground state string (this mode corresponds to the 
Goldstone mode associated to spontaneous symmetry breaking in the Gross-Pitaevskii 
framework), or the two-particle states obtained when splitting it up by extracting a 
single atom. The response of the system to an external probe therefore essentially 
occurs with zero recoil energy, which is interpreted as a Mossbauer-like effect: the 
attractive gas is really like a crystal along the imaginary axis in momentum space, 
and responds like a single particle since the gap to the lowest excitations is finite. For 
large atom number TV, higher families of excited states have their contribution to the 
correlation functions suppressed by progressively higher powers of 1/N, and can be 
ignored for all practical purposes. Our results thus provide an extremely accurate, 
essentially exact picture of the zero-temperature dynamics of the system. 

Using our results, it is possible to compute other more elaborate form factors, for 
example those of the density or field operator between higher excited states involving 
nontrivial string partitioning. These take the form of reduced determinants involving 
only the string centers, but are still obtainable using a reasoning similar to the one 
used in the calculation of the single-string form factors in Appendix C and Appendix 
D. We have not pursued this line here since these would contribute sub-dominantly 
to the observable zero-temperature lineshapes. 

Integrable models for atomic gases also include the class of multi-species, 5- 
interaction particles. For these, the Bethe Ansatz involves a nesting, making the 
current approach somewhat more complicated. It would be interesting to see whether 
the methods used here could be adapted to such cases, and used to derive related 
expressions for correlation functions like the spin structure factor. Such an object 
would also be in principle observable experimentally [58]. 

Another interesting possibility is to consider the effects of a finite temperature. 
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This would probably be of great importance for experiments, since the ground state 
is part of a dispersionlcss band in the large A'^ Umit; any temperature greater than 
0{1/N) would thus quickly smear out the correlations in the momentum direction. 
However, all the necessary form factors for the calculation of the finite-temperature 
response functions are readily obtained directly from the results wc presented, and it is 
more or less straightforward to obtain the finite temperature dynamical correlations. 
The presence of a harmonic trap can also be taken into account by convolving the 
ground state wavepacket with a harmonic oscillator wavefunction. We will investigate 
these issues further in the near future. 

In conclusion, we have shown that the attractive Lieb-Liniger gas is a system for 
which the Bethe Ansatz is sufficiently tractable to allow an analytical calculation of 
important correlation functions. 
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Appendix A. Reduced Bethe equations 

In this appendix we derive the reduced Bethe equations for the string centers. The 
Bethe equation (12) in terms of the string rapidities (15) are 

\j,o. _ \k,b _ ■- 
^iX'^^L ^ TT ^" 

ii \j,a _ , fc,fc , ■- 

(fe,/3,6)#(i,a,a) ^" + ^'^ 

-pr jT K 4 + -a + b-1) c{-a + &-!) + 

iJka) L\ - + ^a + b+1) l[ -ci-a + b+l) + ^^i^-") 

where we have separated inter- and intra-string parts, dropped all string deviations for 
inter-string factors, and denoted S{-" — d{^'' — (5^'^"'^^ in the intra-string part. Simplified 
Bethe equations are obtained by taking the product of these equations within the string 
considered. The left-hand side becomes 
j 

Yl e'^i'"-^ = e'^^i-^. (A.2) 

o=l 

On the right-hand side, we have 

a=ib^ac{-a + b+l) + 6ii' 
for the intra-string part, and (writing A = — A^) 
\ + ic{^- a + b-1) 



^]lx + z-ci^-a + b+l) 
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= e|,_fe|(A)e§_fcj+2(A)e§_fc|+4(A)...e2^fc_2(A)ej+fe(A) = E^kW (A.4) 

where 

X + icj/2 

The exponential form of the Bethe equations has thus been reduced to the set of Ng 
coupled equations for the string centers A^, 

^ijXiL^ n EM-X'ff). (A.6) 

(fe,/3)#(j,c,) 

Taking the logarithm and defining 

2A 

^!>,(A) = 2 atan — (A. 7) 

cj 

such that — i log(— ej(A)) = 4>j{X), we find equations (17) and (18). 



Appendix B. Norm of a general state 



In this appendix, we consider the fate of the Gaudin-Korepin norm formula in the 
case of general string states. A similar treatment can be found in [55]. First of all, 
the prefactor of equation (32) explicitly reads (in string notation) 



n 



{k,(3,b)^{3,a,a) " ^0 



yr A A Aj,-Afe+»c(i^-a + 6+l) 
J=\L\ Ai-A^ + zc-(i-^-a + 6) 



(fe,/3)#(j 



where we have separated the product into inter- and intra-string parts and defined 

_ s:j,a _ ^j,b 
Oa — Oa 0„ . 

To leading order in string deviations, this prefactor is then simplified to (including 

the factor) 

e^3-Q,A^., -Q F,fe(Ai-A^)xnn'^i^"'''^'^x(l + 0('^)) (^-2) 

j (k,l3)>{j,a) j,a o=l 

where 



A^ + (i(j+fc))' 



(B.3) 



Let us now deal with the Gaudin matrix. In string notation, the matrix elements 



read 

Q(j,a,a)(k,0,b) = ^(j,a,a)(k,0,b) 

where 



(j,a,a)(fe,/3,6) 



i+ X! -^(j,a,a)(;,7,c) - K(^j^ct,a){k,l3,b) 
(.1,1,0) 

2c 

■(Ar-A^''')2 + c2' 



(B.4) 



(B.5) 
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The diagonal block associated to a particular string within the Gaudin matrix is 

/ L + K12 + S1 -K12 -Ki3 - -Ki, \ 

K12 L + K12 + K23 + S2 -K23 ■■■ -K2j I 

Ki3 -K23 L + K23 + K34. + S3 ... ... (B.6) 

\ -Kxj ... ... -Kj-ij L + Kj-ij+Sj I 

in which we have isolated elements of the type Kaa+i which are exponentially large, 

All other elements, including the sums Sj (of finite Kab terms) are of zeroth order in 
string deviations. 

The determinant can now be manipulated by adding block column 1 to block 
column 2, followed by adding block row 1 to block row 2. This eliminates the K12 
factor except in the block (1, 1) position. Factors Kaa+i for a, = 2, j — 1 can be 
similarly eliminated from two off-diagonal and one diagonal positions. The string block 
then has j — 1 exponentially large factors Kaa+i on the diagonal, and all off-diagonal 
terms within the block can be neglected. Doing similar operations for all strings in 
the eigenstate and keeping all remaining factors into account, we find that, to leading 
order in the string deviations, the determinant of the Gaudin matrix becomes 

= n n Det^e-) X (1 + 0{6)) (B.8) 

i,a 0=1 "a 

where the reduced Gaudin matrix is an Ng x Ng matrix whose elements are 

,(r) d 



jXiL-Y,^ji{Xi-X\) . (B.9) 

(1,-1) 

Upon taking the product of this with the prefactor (B.2), the string deviation 
products cancel, and we find that the norm of a completely general string state 
(dropping corrections of order 5) is finite and can be written as 

ll*iv({A^c«})ir = c-^=n-?'''^ n Detg^^\ (B.IO) 

j (k,0)>(j,a) 

(r) 

In the limit of large L, we can simplify the reduced Gaudin matrix to G(^j [^^ f^f. p-^ = 
jLSjkSa0, allowing to simplify the norm to (35). 

Appendix C. Density operator form factor 

The form factor of the density operator between two Bethe eigenstates of the Bose 
gas is given in [50, 51] as 

^^A = (P. - P^) -VPII (^^) + U(A,)]^;^^l^ , (C.l) 

with 
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and 



3 3 m^j 



n 



Mr7 



A, 



Am ~ Ai 



(if(A,-Afe)-i^(Ap-Afc)),(C.3) 



and where the kernel K is given by Eq. (34) and Xp an arbitrary number, not 
necessarily a rapidity. 

Fixing Xp = Ajv, the last column {k = N) becomes (0, ...,0, 1) and the form 
factor can be written in terms of an (A'' — 1) determinant: 



N- 



Xjk + ic 



3,k 



Xi 



BetN-iU{{fi},{X}) 



(C.4) 



with the matrix U having elements 

Ujkm, {A}) = Sj /'''~/' + ^l^^^^(K{Xj - A,) - KiX^ - X,)). (C.5) 

Let us now consider the case where rapidities {A} form an A^-string, as for the 
ground state. We leave the set of rapidities {/x} arbitrary. We thus take a simple 
notation in which 



ic , 



Ai = A+-(A^ + l-2i) + i^,-. 



(C.6) 



We have to deal explicitly with the small deviations in our matrices. Namely, we 
have 

^ +0(1), fc = j + l. 



K{X, - Afc) 



1 



+ 0(1), fc=i-l. 



(C.7) 



0(1), 



fc^jil. 



In these, we have written 6 



jk 



5k and kept only the leading term, which turns 



out to be the only one needed for the calculation. 
We also define the regular function 



Ylm=l(^^m - Xj 



0,1,...,7V+1 



(C.8) 



^^jV'rn - Xj) 

by defining Aq and Aat+i using equation (C.6) for j = 0,N+1 (we don't need to define 
6j for these values of j). The Vj^ functions can thus be written in regularized form 



= < 



= < 



Vo, 


3 = 


1, 




i = 2, 


..,N 




J = l, 


..,7V 


Vn+1, 


3 


= N. 



(C.9) 
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To leading order, the U matrix (C.5) is 



V 



_V2 

^12 
V2 








-Vi _ V3 



V3_ 
■523 





V2 



-V2 
523 



V4 
■534 



23. 
534 



-Vl 

"n-in 

-V2 
-V3 



-Vjv_ 



VN-2 



Vn-2 



VN-1 
<>]V-2N-1 



''N-2iV-l 



<'N-2N-1 TJv-lJV 
-VW-2 _ VW-l+ViV 
<>]V-2N-1 ''N'-IN 



(C.IO) 



The determinant can be evahiated by means of simple row and cohimn manipulations. 
Adding the first row to the second one, thereafter the second to the third, as so on, 
we end up with an upper triangular matrix whose determinant is simply 



Det U = (-1) 



N-l 



1 



N-l 

3 %+i 



N-2 

i=2 



N 



(C.ll) 



Using Um^Mm - A,) = (zc)^-i n„^,(-m + j) = {tcr-H-ir-^{j - mN - jy. 

for the denominator of Vj and defining the function 



Qj — J^C/^m ~ Aj 



(C.12) 



for the numerator, we arrive after simple manipulations to the following expression 
for the density operator form factor (taking the ground state momentum to zero) 



and where we have defined the summation 

N 



" -SN,N{{2iHm/c}) (C.13) 



JV 



N- 1 

i-1 



^ M 

J H [N + l-2j+fim]{G.U) 



Sn,m{{iAm) = 

for any set of complex numbers {/v}. This sum can be performed explicitly (see 
Appendix E), yielding the final expression for the density operator form factor between 
the ground state (A''-string centered on zero) and the arbitrary string state defined by 
the set of string rapidities {/U^}: 

NT\N)P^Jc 



yP 



(C.15) 



Appendix D. Field operator form factor 

The form factor of the field operator between two Bethe states is [52, 14] 

T-iN > . 

= ^''-'c-'/' n " , Detiv-it/({A^}, {A}) , (D.l) 

Uj rife Mfe-Aj 

with the matrix U once again given by (C.5), but with the difference that the set of 
rapidities /x now comprises TV — 1 elements. The calculation is almost identical to that 
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for the density operator, except that here we have Qm = TlJ^iil^'j — ^m)^ i-e- it is the 
product of A'' — 1 and not N factors. The prefactor is shghtly different, and we get 

The summation can also be performed explicitly for this new definition of Qj using 
the general results of Appendix E, we finally find 



Note once again that the set of {/z} rapidities comprises — 1 elements instead of N 
like in (C.15), so although these terms look the same, the different jj. string contents 
involved make them different. 

Appendix E. Proof of summation formula 

We wish to compute 

N / W 1 \ ^ 

5;v,m({m}m) = E(-1)'"M ill n[^ + l-2i + Mj (E.l) 

j=l ^ ^ m=l 

for M < and an arbitrary set of parameters {/x}. The product can be written as 

M M 

l[[N + l-2j + = + 1 - 2i)^-'PK{M}) (E.2) 

m=l /=0 

where Pi is the l-th order completely symmetric polynomial. We therefore need to 
compute all 

Tr,,L = J2(-^y~'( ^jll )(iV + l-2i)^ L = 0,...,N. (E.3) 

Consider 

N-l , N 

U^{a) = Y.[-iy ( ) e''(^-i-2^') = e«(^-i)(l-e-2«)^-i = (2sinha)~-i.(E.4) 

This is a generating function for Tjv,z, since 

Tn,l = ^UN{a)\a=o. (E.5) 

For the summations we need in order to calculate the density and field operator form 
factors, we can keep only the leading order in a and write f7jv(a) = 2^~^a^~^ + .... 
Therefore, we have 

Tjv,L = ^L,Jv-i2^-'(iV-l)! (E.6) 

which yields 

Sn,m = 2''-'{N - 1)! (SM^NPiibA) + Sm,n-i) ■ (E.7) 

This can now be directly used in the form factor calculations of Appendix C and 
Appendix D. 



Dynamics of the attractive ID Bose gas: analytical treatment from integrability 25 



References 

A. Gorlitz, J. M. Vogels, A. E. Leanhardt, C. Raman, T. L. Gustavson, J. R. Abo-Shaeer, A. P. 
Chikkatur, S. Gupta, S. Inouye, T. Rosenband and W. Ketterle, Phys. Rev. Lett. 87, 130402 

(2001). 

M. Greiner, I. Bloch, O. Mandel, T. W. Hansch, and T. Esslinger, Phys. Rev. Lett. 87, 160405 
(2001). 

H. Moritz, T. Stofcrlc, M. Kohl and T. EssUngcr, Phys. Rev. Lett. 91, 250402 (2003). 

B. Paredcs, A. Widera, V. Murg, O. Mandel, S. Foiling, I. Cirac, G. V. Shlyapnikov, T. W. 
Hansch and L Bloch, Nature 429, 277 (2004). 

T. Kinoshita, T. Wcnger and D. S. Weiss, Science 305, 1125 (2004). 
T. Kinoshita, T. Wenger and D. S. Weiss, Phys. Rev. Lett. 95, 190406 (2005). 
T. Kinoshita, T. Wenger and D. S. Weiss, Nature 440, 900 (2006). 

I. Bloch, J. Dalibard and W. Zwcrgcr, preprint arXiv;0704.3011. 

E. H. Lieb and W. Liniger, Phys. Rev. 130, 1605 (1963); E. H. Lieb, Phys. Rev. 130, 1616 
(1963). 

T. Giamaxchi, "Quantum Physics in One Dimension" , Oxford University Press, Oxford, 2004. 

F. D. M. Haldane, Phys. Rev. Lett. 47, 1840 (1981). 

A. O. Gogolin, A. A. Nersesyan and A. M. Tsvelik, "Bosonization and Strongly Correlated 
Systems", Cambridge University Press, Cambridge, 1998. 

J.-S. Caux and P. Calabrese, Phys. Rev. A 74, 031605 (2006). 
J.-S. Caux, P. Calabrese, and N. A. Slavnov, J. Stat. Mech. P01008 (2007). 
V. E. Korepin, N. M. Bogoliubov and A. G. Izergin, "Quantum Inverse Scattering Method and 

Correlation Functions" , Cambridge, 1993, and references therein. 
P. Calabrese and J.-S. Caux, Phys. Rev. Lett. 98, 150403 (2007). 
J. B. McGuire, J. Math. Phys. 5, 622 (1964). 

F. Calogcro and A. Dc Gaspcris, Phys. Rev. A 11, 265 (1975). 

Y. Castin and C. Herzog, Comptes Rendus de I'Academie des Sciences de Paris, tome 2, serie 

IV, pp. 419-443 (2001). 
L. Salasnich, J. Phys. B 39, 1743 (2006). 

S. Sinha, A. Yu. Cherny, D. Kovrizhin and J. Brand, Phys. Rev. Lett. 96, 030406 (2006). 

C. Lee and J. Brand, Europhys. Lett. 73, 321 (2006). 

G. M. Kavoulakis, Phys. Rev. A 67, 11601 (2003). 
R. Kanamoto, H. Saito and M. Ueda, Phys. Rev. A 67, 13608 (2003). 

K. E. Strecker, G. B. Partridge, A. G. Truscott and R. G. Hulet, Nature 417, 150 (2002). 
L. Khaykovich, F. Schreck, G. Ferrari, T. Bourdel, J. CubizoUes, L. D. Carr, Y. Castin and C. 
Salomon, Science 296, 1290 (2002). 

B. Eiermann, Th. Anker, M. Albiez, M. Taglieber, P. Treutlein, K.-P. Maxzlin and M. K. 
Oberthaler, Phys. Rev. Lett. 92, 230401 (2004). 

S. Inouye M. R. Andrews, J. Stenger, H.-J. Miesner, D. M. Stamper-Kurn, W. Ketterle, Nature 

392, 151 (1998). 
M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 

T. Bergcman, M. G. Moore and M. Olshanii, Phys. Rev. Lett. 91, 163201 (2003). 
P. D. Drummond, R. M. Shelby, S. R. Friberg and Y. Yamamoto, Nature 365, 307 (1993). 
L.-M. Duan, Phys. Rev. Lett. 96, 103201 (2006). 

J. Stenger, S. Inouye, A. P. Chikkatur, D. M. Stamper-Kurn, D. E. Pritchard and W. Ketterle, 

Phys. Rev. Lett. 82, 4569 (1999). 
T. Stofcrlc, H. Moritz, C. Schori, M. Kohl and T. Esslinger, Phys. Rev. Lett. 92, 130403 (2004). 
V. Gritsev, E. Altman, E. Demler and A. Polkovnikov, Nature Physics 2, 705 (2006). 
A. Polkovnikov, E. Altman and E. Demler, Proc. Nat. Acad. Sci. 103, 6125 (2006). 
Y. Japha, S. Choi, K. Burnett, and Y. B. Band, Phys. Rev. Lett. 82, 1079 (1999). 

D. L. Luxat and A. Griffin, Phys. Rev. A 65, 043618 (2003). 
P. Blair Blackie, New. J. Phys. 8, 157 (2006). 

I. E. Mazets, G. Kurizki, N. Katz, and N. Davidson, Phys. Rev. Lett. 94, 190403 (2005). 
T.-L. Dao, A. Georges, J. Dahbard, C. Salomon, and I. Carusotto Phys. Rev. Lett. 98, 240402 
(2007). 

W. Ketterle, D.S. Durfee, D.M. Stamper-Kurn, in "Proceedings of the 1998 Enrico Fermi summer 
school on Bose-Einstein condensation in Varenna", Italy [cond-mat/9904034]. 

C. N. Yang and C. P. Yang, J. Math. Phys. 10, 1115 (1969). 

H. B. Thacker, Rev. Mod. Phys. 53, 253 (1981). 

M. Takahashi, Thermodynamics of one-dimensional solvable models (Cambridge University 



Dynamics of the attractive ID Bose gas: analytical treatment from integrability 26 



Press, Cambridge, 1999). 

[46] K. Sakmann, A. I. Streltsov, O. E. Alon, L. S. Cederbaum, Phys. Rev. A 72, 033613 (2005). 

[47] A. G. Sykes, P. D. Drummond and M. J. Davis, arXiv:0707.2422. 

[48] M. Gaudin, "La fonction d'onde de Bethe", Masson, Paris, 1983. 

[49] V. E. Korepin, Commun. Math. Phys. 86, 391 (1982). 

[50] N. A. Slavnov, Teor. Mat. Fiz. 79, 232 (1989). 

[51] N. A. Slavnov, Tcor. Mat. Fiz. 82, 389 (1990). 

[52] T. Kojima, V.E. Korepin, and N.A. Slavnov, Commun. Math. Phys. 188, 657 (1997). 

[53] C. Astrakharchik, J. Boronat, J. CasuUcras, and S. Giorgini, Phys. Rev. Lett. 95, 190407 (2005). 

[54] M. T. Batchclor, M. Bortz, X. W. Guan, and N. Oclkcrs, ,1. Stat. Mech. LlOOOl (2005). 

[55] A. N. Kirillov and V. E. Korepin, J. Math. Sci. 40, 13 (1988). 

[56] M. Olshanii and V. Dunjko, Phys. Rev. Lett. 91, 090401 (2003). 

[57] A. Minguzzi, P. Vignolo and M. P. Tosi, Phys. Lett. A 294, 222 (2002). 

[58] L Carusotto, J. Phys. B 39, 8211 (2006). 



